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Abstract 

We place observational constraints on a coupling between dark energy and dark matter by using 71 Type la 
supernovae (SNe la) from the first year of the five-year Supernova Legacy Survey (SNLS), the cosmic microwave 
background (CMB) shift parameter from the three- year Wilkinson Microwave Anisotropy Probe (WMAP), and the 
baryon acoustic oscillation (BAO) peak found in the Sloan Digital Sky Survey (SDSS). The interactions we study 
are (i) constant coupling 5 and (ii) varying coupling 5{z) that depends on a redshift z, both of which have simple 
parametrizations of the Hubble parameter to confront with observational data. We find that the combination of 
the three databases marginalized over a present dark energy density gives stringent constraints on the coupling, 
-0.08 < 5 < 0.03 (95% CL) in the constant couphng model and -0.4 < (5o < 0.1 (95% CL) in the varying couphng 
model, where 5q is a present value. The uncoupled ACDM model [wx = — 1 and (5 = 0) still remains a good fit to the 
data, but the negative coupling {5 < 0) with the equation of state of dark energy wx < — 1 is slightly favoured over 
the ACDM model. 



PACS number(s): 98.80.Es, 98.80. Cq 

I. INTRODUCTION 

Recent observations of Supernova la (SNIa) suggest 
that the universe has entered a stage of an accelerated 
expansion with a redshift z ^1 [H, 0j 13] ■ This has been 
confirmed by precise measurements of the spectrum of 
the Cosmic Microwave Background (CMB) anisotropies 
[1, [H] as well as the baryon acoustic oscillations (BAO) 
in the Sloan Digital Sky Survey (SDSS) luminous galaxy 
sample @. As is well known, all usual types of matter 
with positive pressure generate attractive forces, which 
decelerate the expansion of the universe. Given this, a 
dark energy component with negative pressure was sug- 
gested to account for the invisible fuel that drives the cur- 
rent accelerated expansion (see Refs. for reviews). 

The simplest candidate for dark energy is a cosmolog- 
ical constant A (vacuum energy), which corresponds to a 
constant equation of state wx = — 1- This model, the so- 
called ACDM, provides an excellent fit to a wide range of 
astronomical data so far. However, in such a model there 
exists a theoretical problem of cosmic coincidence: why is 
the vacuum density comparable with the critical density 
at the present epoch in the long history of the universe? 
One possible approach to alleviating this problem is to 
assume that the "cosmological constant" is not a constant 
but is a dynamical component with a slowly evolving and 
spatially homogeneous scalar field called quintessence [9] 
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(see Refs. [T^ for early works and Refs. [Tl| for the re- 
construction of quintessence potentials) . In such models 
the resolution of the cosmic coincidence problem typically 
leads to a fine-tuning of model parameters. 

Given the fact that the amount of dark matter is 
comparable to that of dark energy in the present uni- 
verse, it is natural to consider an interaction between the 
two components. Originally cosmological consequences 
of a scalar field coupled to the matter were studied in 
Ref. [13 • Amendola [3] considered interaction between 
a quintessence field </> and dark matter with a coupling 
Q that satisfies the relation T'^^^y^ = QT(^^^(j).^^ between 
energy momentum tensors. In fact this type of inter- 
action appears in the context of scalar-tensor theories 
[l^ . f{R) gravity model s [l5| . varying mass dark mat- 
ter/neutrino models p^ . Il7l. Il8l. ITgj and phantom dark 
energy models [2^. Interestingly the equation of state 
of dark energy can extend to the region wx < — 1 if 
the mass of dark matter depends on a quintessence field 
[itL [TsI I . The presence of the coupling can also provide 
an accelerated scaling attractor along which the ratio of 
the energy densities of dark energy and dark matter is 
a constant. Then this may be useful to solve the co- 
incident problem because the present universe can be a 
global attractor with a dark energy fraction fix — 0.7. 

However, the coupling Q to realize the accelerated scal- 
ing attractor is too large to ensure the presence of a 



2 



standard matter dominated epoch In fact there ex- 
ists the so-called "(/> matter-dominated epoch" ((/)MDE) 
during which an effective equation of state is given by 
w^g = 2QV3 B E Amendola ^ showed that 

the coupled quintessence model with an exponential po- 
tential is not consistent with observational data of CMB 
unless the coupling Q is smaller than the order of 0.1, 
but in this case there is no scaling accelerated attractor. 
Instead the system finally approaches a scalar-field dom- 
inated attractor with Wcs = — 1 and f2x = 1, in which 
case the coincident problem is not solved. 

There are many other scalar-field dark energy models 
such as K-essence, tachyon, phantom and dilatonic ghost 
condensates. For a general Lagrangian density X) 
with a kinetic term X = — (1/2)(V0)^ and a constant 
coupling Q, the existence of scaling solutions required 
to solve the coincident problem restricts the Lagrangian 
density to the form p — X g{Xe^'^), where g is an ar- 
bitrary function and A is a constant [23 |. Recently it 
has been shown that for the vast class of this generalized 
Lagrangian that include most of scalar-field dark energy 
models, the matter era is not followed by the accelerated 
scaling attractor [1^. Thus it is still a challenging task 
to construct coupled scalar-field models which can solve 
the coincidence problem without using fine-tuned varying 
coupfings [22|. 

Since the origin of dark energ y is not yet known, there 
are several different approaches [251 . [26l . [27l . [28', 29] to im- 
plementing couplings without restricting to scalar-field 
models (see also Refs. [33|-[5^ for a number of interest- 
ing aspects of interacting dark energy). The approach 
we adopt in this paper is to introduce an interaction of 
the form Tpm on the rhs of conservation equations, see 
(21) and ([3]). While this is basically a fiuid description 
of dark energy, Eqs. ([2]) and ([3]) include the aforemen- 
tioned scalar-field coupling by setting F — Qcj). Since 
the interaction rate F measured by the Hubble rate H 
is generally important to discuss the strength of an en- 
ergy transfer, we introduce a dimensionless coupling 6 in 
the form S = T/H. This is an approach a number of 
authors adopted [1^. One can constrain the strength of 
the interaction observationally by assuming that 5 is a 
constant (as in the constant Q case discussed above). In 
fact the authors in Ref. f29(| recently placed observational 
constraints on the coupling by using the SNIa data with 
a parametrization of the dark energy equation of state 

wx = Wo + WiZ. 

It is also possible to address the varying 5 case. For ex- 
ample, Dalai et al. [2^ assumed that the ratio of dark en- 
ergy and dark matter has a relation px / Pm oc , where 
a is a scale factor and ^ is a constant. For a constant 
equation of state of dark energy, the coupling 5 is known 
as a function of the redshift z [23, see Eq. (fT2)) . Since 
the strength of the coupling decreases for larger z, this 
scenario can address the situation in which the interac- 
tion is weak during the matter era but becomes strong 
in the dark energy dominated epoch. 

Observational constraints on the coupling 6 — T/H 



have been obtained by using SNIa data [13, ■ In this 
paper, we carry out likelihood analysis of coupled dark 
energy models by using 71 high-redshift SNe la from the 
first year of the five-year SNLS, the CMB shift parameter 
from the three-year WMAP observations and the BAO 
peak found in the SDSS. We concentrate on two classes 
of interacting models: (i) a constant coupling S and (ii) 
a varying coupling S{z) with the relation px/pm oc . 
Throughout this paper, the dark energy equation of state 
Wx is assumed to be a constant. In both models, we 
have three free parameters {S^wxt^xo), where l^xo is 
the present energy fraction of dark energy (in the varying 
coupling model 5 is replaced by the present value 5q). We 
find that the combination of the three databases places 
stringent constraints on the model parameters since the 
CMB shift parameter is sensitive to the value of the cou- 
pling. Our results also indicate that the concordance 
ACDM model still remains a good fit to the data, but 
the negative coupling ((5 < 0) with wx < — 1 is slightly 
favoured over the ACDM model. 



II. INTERACTIONS BETWEEN DARK 
ENERGY AND DARK MATTER 

In this section, we explain the form of the interaction 
between dark energy and dark matter. The background 
metric is described by the flat Friedmann-Robertson- 
Walker (FRW) metric with a scale factor a: 

ds^ = -d<2 + a2(i)dx2, (1) 

where i is a cosmic time. Quite generally we can write 
the conservation equations in the forms 

Pm + 3i//9„i = +Fp„i , (2) 
Px + 'iH{px +px) = -'^Pm , (3) 

where H = a/a is a Hubble rate, pm and px are the 
energy densities of dark matter and dark energy respec- 
tively, and Px is the dark energy pressure density with 
the equation of state wx = Px/px- If we consider a 
scalar-field model of dark energy, the interaction term is 
typically given by F = Q^, where the constant Q char- 
acterizes the strength of the coupling [l^, [2l| . Amendola 
obtained the constraint on the coupling as Q < 0.08 by 
using the information of the CMB power spectrum. 

We note that the origin of dark energy is not yet iden- 
tified as a scalar field. In this work we take a different 
approach to constraining the strength of the interaction 
without assuming scalar-field models. We measure F in 
terms of the Hubble parameter H and define the dimen- 
sionless coupling 

5 = T/H. (4) 

Note that a positive 5 implies a transfer of energy from 
dark energy to dark matter, and vice versa. From 
Eqs. (HI) and ([3]) it is clear that the total energy density is 
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conserved. Neglecting the contributions of (uncoupled) 
baryon and radiation components the Fricdmann equa- 
tion is given by 

3i/2 = K2(p,„ + px), (5) 

where — 8ttG with G being gravitational constant. 

In the following, we first discuss the case in which 6 is 
a constant. Then the analysis is extended to the case in 
which S varies in time and the ratio of the energy densities 
of dark energy and dark matter scales as px / Pm oc • 



A. Constant coupling models 



For constant 6 [H, [2^ , Eq. ^ is easily integrated to 




Pm — Pmoa 



-3+5 



p„io(l + z) 



3-5 



(6) 



where the subscript "0" represents the present values. 
Note that z is a redshift which is defined by z = ap/a — 1 , 
where oq is normalized as oq = 1. Equation ([6]) shows 
that the interaction leads to a deviation from the usual 
conservation relation pm oc a~^. 

We assume that wx is a constant. Then substituting 
Eq. ([6]) into Eq. ([3]), we obtain the following integrated 
solution 



FIG. 1: Probability contours from the SNLS data only at 
68.3%, 95.4% and 99.7% confidence levels in the {wx, S) plane 
marginalized over Qxo with priors ^Ixo ~ 0.72 ± 0.04 and 
(5 < 3 in the constant coupling models. The horizontal and 
dashed lines represent the uncoupled "XCDM" models and 
the coupled ACDM models respectively, and their crossing 
point corresponds to the standard ACDM model. In this case 
we have the constraint —1.78 < 5 < 3 (95% CL). 



Px = pxoil + 

S 



-PmO 



6 + 3wx 



(7) 



Using the Friedmann equation ([5]), we find 



E\z) = nxo{l + z) 



3{l+wx) 



1 - Q 



xo 



S + Swx 



(5(1 + ^)3(1+"'^) +3u;x(l + z) 



3-S 



(8) 



where E{z) = H{z)/Hq and f^xo = K^pxo/{iH^)- Thus 
we have three free parameters (5, wx, ^xo) when we con- 
front models with observations. This allows us to param- 
eterize a wide range of possible cosmologies in a simple 
fashion. 

In the high redshift region {z » 1), it follows from 
Eq. d?]) that px behaves as px — —PmO <5/ {S + 3wx){^ + 
z)^^^ for 3wx < —S. This means that the energy density 
of dark energy becomes negative for S < 0. Since such 
a negative energy appears in phantom models [S^l and 
also modified gravity models [3, UHl , we do not exclude 
the possibility of the negative coupling. Note also that 
in the high redshift region we have the scaling relation 
Px/ Pm = —5/ {5 + iwx)- If \5\ is much smaller than 1, 
the ratio satisfies the relation \px / Pm \ ^ 1 provided that 
Wx is of order —1. 

With the parametrization given above, we can classify 
the models into the following four types in the {wxi5) 



plane: (i) decaying phantom characterized by (5 > and 
Wx < —1, (ii) decaying quintessence"* characterized by 
S > and wx > ^ 1, (iii) created quintessence character- 
ized by (5 < and wx > —1, and (iv) created phantom 
characterized by (5 < and wx < —1, as shown in Fig. [T] 
In this plane, the horizontal dashed line represents the 
uncoupled models with 6 = and the vertical dashed line 
represents the coupled ACDM models with wx = — 1. 



B. Varying coupling models 

When S varies in time, Eq. ^ together with Eq. ^ 
gives 



0-0 na-^eZ-^dOna) 



(9) 



Let us now consider a situation in which the ratio of dark 
energy and dark matter has the following relation [25] : 



Px _ Pxo 5 

Pm PmO 



(10) 



where ^ is a constant which quantifies the severity of 
the coincidence problem. In the absence of the coupling 
6 with constant wx, the energy density of dark energy 



* Here we use the word "quintessence" for dark energy satisfying 
Wx > — 1 without restricting to scalar-field models. 
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scales as px oc a~'^^^~^'^^\ Here the ratio px/pm is 
proportional to a~'^™^, namely, the ^ = — 3wx case in 
Eq. (Uni). Note that the standard ACDM model corre- 
sponds to wx — —1 and £, — '?>■ 

The general case ^ 7^ — Suijis: indicates the existence 
of an interaction between dark matter and dark en- 
ergy. Using the relation (fTUl) . we find that the cou- 
pling r is given by F = —H{£, + 3wx)^xiz), where 
^x{z) = Px / {iH'^)- This shows that 5 varies accord- 
ing to the change of Q.x as 



5{z)^-{i + ?>wx)^x{z). 



(11) 



When ^ = —iwx, this reduces to (5 = 0. Since Vlx{z) is 

given by VLx{z) = [(p,„o//Oxo)(l + z)^ + l] under the 
condition PU)) . the coupling can be written as 



So 



flxo + {I - ^xo){l + z)i 



(12) 



where 60 = — -I- 3wx)^xo- If C > 0; '^(^) decreases for 
higher z. 

Note that the ^ = case gives a constant coupling 
5{z) = Jo- Since this corresponds to an exact scaling 
solution px oc pm-, one cannot realize the matter dom- 
inated epoch followed by a late-time acceleration. This 
constant 6 case is different from the one we discussed in 
subsection A. In fact the solutions ([6]) and ([7]) do not 
satisfy the relation (fTU)) . 

From Eqs. © and © together with Eq. the total 
energy density pr — Pm + Px satisfies 



d In pt 
da 



1 I f -e 

1 + wx a ^ 

\Pxo 



which can be integrated to give 



Pt = Pro a 



Pxo 
Pto 



(l-a«) 



(13) 



(14) 



Then the Friedmann equation ([5]) gives 



E^z) = (1 + zf [I - nxo + ^xoil + 

(15) 

From Eq. (fTO|) . we find that the energy density of dark 
energy is given by 



Px = Pt 



^xa 



rtxo + {I - i^xo){l + z)i 



(16) 



We have three parameters (^, wx, ^xo) in this model. 
Since 60 is related to these variables by the relation 5q = 
— (^ -|- 3wx)^X0j one can instead vary three parameters 
{Sq, Wx, ^xo) when we carry out likelihood analysis. 
Thus the model (jlSp is a simple parametrization that 
implements the variation of the coupling S. Note that 
both Pt and px are positive as long as < flxo < 1- 
Hence px remains positive in the high-redshift region 
even for 6 < unlike the constant coupling model. 



The coupling affects the 
of interest, such as the age 
eration parameter. Given 
universe and the transition 
switches from deceleration 
as the value of Sq increases 
the next section, we place 
the strength of the couplinj 



evolution of some quantities 
of the universe and the decel- 
wx and flxo, the age of the 
redshift at which the universe 
to acceleration become larger 
from negative to positive. In 
observational constraints on 



III. CONSTRAINTS FROM RECENT 
OBSERVATIONS 

In this section, we study the viability of interacting 
models presented in the previous section by using re- 
cently released SNLS data Q in conjunction with the 
BAG peak in the SDSS [1] and the CMB shift parame- 
ter [561. 

Recently Astier et al. Q have compiled a new sample 
of 71 high-redshift SNe la, in the redshift range 0.2 < 
z < 1.0, discovered during the first year of the 5-year 
SNLS. This data set is arguably the best high-redshift 
SN la compiled data, since the multi-band, rolling search 
technique and careful calibration are adopted. The lumi- 
nosity distance ^^(z) to supernovae is given by 



dLiz) = H^^l + z) 



dz' 

W) 



(17) 



The baryon oscillations in the galaxy power spectrum 
are imprints from acoustic oscillations prior to recombi- 
nation, which are also responsible for the acoustic peaks 
seen in the CMB temperature power spectrum. The 
physical length scale associated with the oscillations is 
set by the sound horizon at recombination, which can 
be estimated from the CMB data Measuring the 
apparent size of the oscillations in a galaxy survey al- 
lows one to measure the angular diameter distance at 
the survey redshift. Although the acoustic features in 
the matter correlations are weak on large scales, Eisen- 
stein et al. 0] have successfully found the peaks using a 
large spectroscopic sample of luminous red galaxies from 
SDSS e5_71]. This sample contains 46,748 galaxies covering 
3816 square degrees out to a redshift of 2; = 0.47. They 
found a parameter A, which is independent of dark en- 
ergy models [6J- From Eq. (5) in their paper (fj], we write 
it as 



A = y/a^E{zi) 



-1/3 



1 

zi Jo 



E{z') 



2/3 



(18) 



where zi — 0.35 and A is measured to be A = 0.469 ± 
0.017. In our analysis, we combine these measurements. 

The CMB shift parameter TZ captures the correspon- 
dence between the angular diameter distance to last scat- 
tering surface and the relation of the angular scale of 
the acoustic peaks to the physical scale of the sound 
horizon Its value is expected to be mostly model- 

independent, which can be extracted accurately from 
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CMB data. The shift parameter TZ is given by [s^ 

where Zj-ec is the redshift of recombination. It provides 
a useful constraint on evolving dark energy models since 
the integral over E{z) extends to high redshifts. The 
recent analysis of the three-year WMAP data [H] gives 
7^ = 1.70 ± 0.03 at Zrec = 1089 

A. Constant coupling models 

Let us first consider observational constraints on con- 
stant coupling models. In Fig. [H we show probability 
contours from SNLS when we take a prior for fJxo such 
that the probability distribution is a Gaussian with a 
mean of 0.72 and a standard deviation given by ct = 0.04. 
We also assume the condition (5 < 3 to exclude the pos- 
sibility that dark matter behaves as a phantom matter 
[see Eq. The phantom dark matter is problematic 

for successful structure formation. Then we find that the 
SNLS data gives a weak constraint on the coupling as 
— 1.78 < 5 < 3 at the 95% confidence level. When we 
choose the wider range of the prior for Qxo, the con- 
straint becomes weaker. 

Figure[2]shows the case in which the BAO data is taken 
into account in addition to the SNLS data without a prior 
for r^xo- Compared to Fig. [U the allowed range of wx is 
reduced. However we still have a rather large region of a 
parameter space for S, i.e., —1.73 < 6 < 3 (95% CL). 

We have also carried out a likelihood analysis without 
any prior for S and found that even the large coupling 
such as 5 = 20 with wx ^ —0.7 is within the 2a obser- 
vational contour bound. This reflects the fact that the 
SNIa and BAO data give constraints only around low 
redshifts z < 0(1). The models can fit the data even for 
i5 ^ 1 because of the dominance of the S{1 + 
term instead of the usual 3wx(l + z)^^^ term on the rhs 
of Eq. ([S]). This tells us how it is important to include 
other data in a high-redshift region (z ^ 1) in order to 
rule out models with problematic couplings {6 > 3). 

In fact the CMB shift parameter provides a stringent 
constraint on the coupling. In Fig. [3] we plot observa- 
tional contours from the joint analysis of SNLS, CMB 
and BAO data in the {wx , S) plane marginalized over 
^xa (left panel) and in the {flxo,S) plane marginalized 
over Wx (right panel). Note that we do not put any 
prior for S in these analysis. We find that the coupling 
is severely constrained: -0.08 < S < 0.03 (95% CL) in 
both marginalizations. This comes from the fact that 
the CMB data do not allow a large deviation from the 
standard matter-dominated epoch. 

The combined analysis of three databases constrains 
the equation of state and the present energy fraction of 
dark energy to —1.16 < wx < —0.91 and 0.69 < flxo < 
0.77 (95% CL). It is interesting to note that the allowed 
observational contours are widely spread in the phantom 
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FIG. 2: Probability contours from the SNLS and BAO data in 
the {wx,S) plane marginalized over ^Ixo without a prior for 
Q,xo and with a prior 5 < 3 in the constant coupling models. 
In this case we have the constraint — 1.73 < 5 < 3 (95% CL). 



region (wx < —1) with a negative coupling (S < 0), 
see Fig. [3] The best-fit parameters are found to he S — 
-0.03, wx = -1.02 and Qxo = 0.73 with = 60.94, 
which is slightly favoured over the ACDM model. 



B. Varying coupling models 

We now proceed to the varying coupling models in 
which S depends upon z in the form (|12p . In Fig. 3] we 
show observational contours from the combined analy- 
sis of SNLS+BAO+CMB data in the planes (i) {wx,So) 
marginalized over flxo (left panel) and (ii) {^IxotSo) 
marginalized over wx (right panel). We find that the 
present coupling, the dark energy equation of state and 
the present dark energy density are constrained to be 
-0.4 < So < 0.1, -1.18 < wx < -0.91 and 0.69 < 
^xo < 0.77 at the 95% confidence level. The best-fit 
parameters correspond to S = —0.11, wx = —1.03 and 
^xo — 0.73 with = 60.94. Similarly to the constant 
coupling case, the SNIa and BAO data do not provide 
stringent constraints on 60, but inclusion of the CMB 
data significantly reduces the allowed region of the cou- 
pling. Since S{z) decreases for larger z, the observational 
constraints on 60 is not so severe compared to the con- 
stant coupling models. We also find that the phantom 
models {wx < — 1) with a negative coupling {60 < 0) 
have a wider allowed parameter space compared to three 
other divided regions in the left panel of Fig. 21 

As we see from Fig. [31 the ACDM model, which corre- 
sponds to the point {wx,S,) = (—1,3), is within the la 
contour bound. We remind that the uncoupled models 
are characterized by the line ^ = —5wx- Thus, provided 
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FIG. 3: Probability contours from the combination of SNLS, BAO and CMB data in the constant couphng models. The left 
panel shows observational contours in the {wx, S) plane marginalized over Qxo without a prior for Qxo, whereas the right panel 
shows contours in the {Qxo,S) plane marginalized over wx with no prior for wx- The best-fit model parameters correspond 
to 5 = -0.03, wx = -1.02 and ^Ixo = 0.73 with = 60.94. In this case we have the constraint -0.08 <5 < 0.03 (95% CL). 



+ BAO + CMB 



0.65 0.7 0.75 0.8 0.85 

FIG. 4; Probability contours from the joint analysis of the SNLS, BAO and CMB data in the varying coupling models. The 
left panel shows observational contours in the (wxjSo) plane marginalized over Qxo without prior for Qxo, whereas the right 
panel shows contours in the {Qxo, S) plane marginalized over wx without prior for wx- The best-fit parameters correspond to 
5 = -0.11, Wx = -1.03 and Qxo = 0.73 with = 60.94. In this case we have the constraint -0.4 < So < 0.1 (95% CL). 



that the points are not on the line ^ = —3wx, the cou- IV. CONCLUSIONS 

pled models are allowed observationally in the parameter 
regions 2.66 < < 4.05 (95% CL). From Fig. [5] it is obvi- 
ous that the scaling models with ^ = are excluded from In this paper, we have studied interacting models of 
the data. dark energy in which the coupling 6 between dark mat- 

ter and dark energy is given by ([4]) together with con- 
servation equations ^ and ^ . We discussed two differ- 
ent models: (i) constant 6 case, and (ii) varying S case 
which has a redshift dependence given in IT^ . The lat- 
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FIG. 5: Probability contours in the varying coupling models 
in the (wx,^) plane marginalized over Qxo- The line ^ = 
—3wx corresponds to the uncoupled models. In this case we 
have the constraint 2.66 < ^ < 4.05 (95% CL). 



ter type of couplings appears by imposing the relation 
Px/Pm oc between the energy densities of two dark 
components. Assuming that the equation of state of dark 
energy wx is a constant, we obtain the convenient forms 
of the Friedmann equations ^ and (|15p to confront with 
observational data. 

We have placed observational constraints on the 
strength of the coupling S by using the recent data of 
SNLS, the CMB shift parameter and the SDSS baryon 
acoustic oscillations (BAO). In both constant and vary- 
ing coupling models, the supernova data alone do not 
provide stringent constraints on the coupling. Adding 
the BAO data generally reduces the allowed range of wx , 
but the coupling S of order unity is not still ruled out. 
This reflects the fact that the BAO and SN la data are 
sensitive to the value of wx or flxo rather than 6. How- 
ever the combined analysis of SNLS-I-BAO-I-CMB shows 
that the allowed region of the coupling is significantly 
reduced compared to the case without the CMB data. 
This is associated with the fact that a large coupling 
leads to the change of the cosmological evolution during 



the matter-dominated epoch, thus modifying the CMB 
angular-diameter distance. 

By the joint analysis of SNLS-I-BAO+CMB, we have 
obtained observational constraints on the strength of the 
couphng: (i) -0.08 <6 < 0.03 (95% CL) for the constant 
couphng models and (ii) -0.4 < (5o < 0.1 (95% CL) for 
the varying coupling models (here So is a present value). 
We also find that the best-fit values exist in the phantom 
region {wx < —1) with a negative coupling {S < 0) in 
both constant and varying coupling models. The uncou- 
pled ACDM model {wx = — 1 and (5 = 0) still remains 
a good fit to the data. Nevertheless it is interesting to 
note that the negative coupling with a dark energy with 
Wx < — 1 is slightly favoured over the ACDM model. 

There exist coupled dark energy models [13, [3 that 
give Wx < — 1 without using a negative kinetic energy 
of a scalar field (called "super-acceleration" in Ref. 18] ). 
This can be realized by considering scalar-field dependent 
masses of dark matter particles. It was shown in Ref. [l^ 
that this super-acceleration model satisfies a number of 
observational constraints, which is consistent with our 
results that the equation of state wx < — 1 is favoured. 

In this work we did not take into account observational 
constraints from matter density perturbations Sm which 
is also affected by the presence of the coupling j^gfl,]. It 
was shown in Ref. [6l| that the inclusion of matter den- 
sity perturbations can place a stronger constraint on the 
coupling S compared to the analysis using the background 
evolution only. The future galaxy surveys such as KAOS 
and PANSTARS will further provide good data of the 
matter power spectrum with a high accuracy, which can 
offer a possibility to put severe constraints on the cou- 
pling. This may provide us an exciting possibility to 
reveal the origin of dark energy and dark matter. 
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